Adiabatic pumping through a quantum dot in the Kondo regime: 
Exact results at the Toulouse limit 



O 

o 

(N 

> 
O 

CO 



Avraham Schiller 1 and Alessandro Silva 2 
1 Racah Institute of Physics, The Hebrew University, Jerusalem 91904, Israel 
2 Abdus Salam ICTP, Strada Costiera 11, 34100 Trieste, Italy 

Transport properties of ultrasmall quantum dots with a single unpaired electron are commonly 
modeled by the nonequilibrium Kondo model, describing the exchange interaction of a spin- 1/2 
local moment with two leads of noninteracting electrons. Remarkably, the model possesses an exact 
solution when tuned to a special manifold in its parameter space known as the Toulouse limit. We use 
the Toulouse limit to exactly calculate the adiabatically pumped spin current in the Kondo regime. 
In the absence of both potential scattering and a voltage bias, the instantaneous charge current is 
strictly zero for a generic Kondo model. However, a nonzero spin current can be pumped through the 
system in the presence of a finite magnetic field, provided the spin couples asymmetrically to the two 
leads. Tunneling through a Kondo impurity thus offers a natural mechanism for generating a pure 
spin current. We show, in particular, that one can devise pumping cycles along which the average 
spin pumped per cycle is closely equal to %. By analogy with Brouwer's formula for noninteracting 
systems with two driven parameters, the pumped spin current is expressed as a geometrical property 
of a scattering matrix. However, the relevant scattering matrix that enters the formulation pertains 
to the Majorana fermions that appear at the Toulouse limit rather than the physical electrons that 
carry the current. These results are obtained by combining the nonequilibrium Keldysh Green 
function technique with a systematic gradient expansion, explicitly exposing the small parameter 
controlling the adiabatic limit. 

PACS numbers: 72.15.Qm, 72.10.Fk 
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I. INTRODUCTION 

The act of pumping is well known from everyday life. 
By repeatedly operating a periodic sequence of steps one 
can transfer a certain amount of fluid or gas between 
reservoirs held at equal potential. The same principle 
applies to electrical charge. By periodically modulating 
spatially-confined potentials it is possible to generate a 
nonzero dc current between leads that are kept at equal 
temperature and electrochemical potential. When oper- 
ated sufficiently slow such that the typical scattering time 
for electrons is much faster than the time over which the 
scattering potentials vary, this process is know as adia- 
batic quantum pumping. Recently, there has been con- 
siderable theoretica^ ^^^^ 12 i 13 i 14 i 15 i 16 i 17 i 18 i 19 
and experimenta l 20 ' 21 ' 22 ' 23 interest in adiabatic quantum 
pumping in confined nanostructures. Besides the fun- 
damental and technological importance of understanding 
time-dependent phenomena in nano-devices such as semi- 
conductor and carbon nanotube quantum dots, adiabatic 
quantum pumping offers new possibilities that otherwise 
are difficult to realize in conventional dc transport mea- 
surements with a finite voltage bias. Most notably, the 
ability to pump a quantized amount of charge per cy- 
cle ; 20 ' 22 which is of potential metrological importance. 
In this paper we address another such example, the gen- 
eration of pure, possibly quantized spin current without 
any charge current ^Uil^MLl^i^ 

In the absence of interactions, adiabatic pumping is 
by now well understood. In particular, building on the 
scattering approach of Biittiker et aL^i Brouwer has ele- 
gantly shown 2 - that the adiabatically pumped current can 



be expressed in terms of the instantaneous (equilibrium) 
scattering matrix. In the case of two driven parameters, 
the pumped charge per cycle reduces to a geometrical 
property of the equilibrium scattering matrix, pertaining 
to the area enclosed in parameter space by the pumping 
cycle. All other details of the pumping cycle, i.e., the ex- 
plicit time dependences of the scattering potentials, are 
irrelevant as long as pumping is adiabatic. 

Far less understood are the effects of interactions, 
where efforts have focused thus far on zero— ' 11 ' 13 and 
one-dimensiona l 12 ' 14 systems. The difficulty with incor- 
porating interactions lies in the need to treat retarda- 
tion effects beyond the static limit. Indeed, recent at- 
tempts to generalize Brouwer's formula so as to include 
interactio n] 25 ' 26 ' 27 have either required the introduction 
of complicated vertex corrections ; 25 ' 26 or the applica- 
tion of a gradient expansion to interaction-induced self- 
energies^ Both formulations can only be implemented 
approximately at this stage, urging the need for bench- 
mark results against which approximate treatments can 
be tested. In this paper we provide such an exact result 
for the pumped currents through a Kondo impurity. 

Kondo-assisted tunneling has been observed by now 
in an abundance of nanostructures, ranging from semi- 
conductor— and nanotube 2 ^ quantum dots, to single- 
atom^ , and single- molecule 3 ^ transistors. In the Kondo 
regime, these systems are described by the well-known 
Kondo model: a spin-i local moment undergoing antifer- 
romagnetic spin-exchange with the conduction electrons 
in the leads. The nonequilibrium Kondo model, either 
with a static or a time-dependent voltage bias, is a dif- 
ficult problem. Remarkably, it possesses an exact solu- 
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tion when tuned to a special manifold in its parameter 
space, known as the Toulouse limit i 32 i 33 At the Toulouse 
limit one can apply a suitable canonical transformation 
to recast the interacting problem in free, quadratic form. 
This requires the introduction of new fermionic degrees 
of freedom having no simple relation to the physical elec- 
trons in the leads. The resulting solution, which gener- 
alizes previous exact results for the equilibrium Kondo 
probleni ) 34 i 35 does not correspond to realistic parame- 
ters. It requires large values of certain exchange cou- 
plings (see below), rendering it incapable of describing 
weak-coupling physics. However, the Toulouse limit is 
expected to correctly describe the strong-coupling regime 
of the nonequilibrium Kondo effect, as different micro- 
scopic models are governed by the same strong-coupling 
fixed point. Indeed, previous applications of the model 
to dc ) 32 i 33 ac,— and pulsed-bias potentials 3 ^ have shown 
all the qualitative features of Kondo-assisted tunneling: 
a zero-bias anomaly that splits in an applied magnetic 
field; Fermi-liquid characteristics in the low-T and low- 
V differential conductance; side peaks in the differential 
conductance at eV = ±ntuv for an ac drive of frequency 
u; and a hierarchy of time scales for the rise, saturation 
and falloff of the current in response to a pulsed bias po- 
tential. The Toulouse limit was also recently applied to 
compute the full counting statistics for tunneling through 
a Kondo impurity. 38 ! 39 

In this paper we take the solution one step further, by 
exactly computing the adiabatically pumped currents on 
the Toulouse manifold. Contrary to previous applications 
of the Toulouse limit to Kondo-assisted tunneling we set 
the voltage bias to zero, but consider a general periodic 
modulation of the transverse exchange couplings and the 
local magnetic field (the free parameters on the Toulouse 
manifold) . In the limit of slow time variations we obtain 
an exact analytic expression for the adiabatically pumped 
spin current. In particular, we show that a nonzero spin 
current can be pumped through the system for a time- 
varying magnetic field, provided the couplings to the two 
leads are made asymmetric. Such a condition is easily 
met in practical devices. Unlike the spin current, the in- 
stantaneous charge current is strictly zero in the absence 
of both potential scattering and a voltage bias, as follows 
from general symmetry considerations. This feature is 
generic to the Kondo model, independent of the adiabatic 
and Toulouse limits. Hence tunneling through a Kondo 
impurity offers a natural mechanism for the realization of 
a spin battery, i.e., a source of pure spin current without 
any charge current. This statement, valid both in the adi- 
abatic limit and beyond, is in qualitative agreement with 
earlier slave-boson mean-field studies of adiabatic pump- 
ing through an Anderson impurity ) 13 ! 25 indicating that 
no fine tuning of model parameters is required as long 
as one operates in the Kondo regime. Finally, we show 
that one can devise suitable pumping cycles that operate 
as a quantized spin pump. Namely, a spin closely quan- 
tized to h is pumped per cycle without an accompanying 
charge. 




FIG. 1: Schematic description of the physical system. A spin- 
| local moment f is placed in between two leads of noninter- 
acting spin-i electrons. The local moment f experiences a 
spin-exchange interaction with the local conduction-electron 
degrees of freedom near the junction, as described by the 
Hamiltonian of Eq. {TJ. Tunneling between the leads is facil- 
itated by spin-exchange terms that scatter an electron across 
the junction. 



As indicated above, the solution at the Toulouse 
limit relies on a nonlocal transformation that converts 
the original spin-exchange Hamiltonian to free-fermion 
form i 33 ' 35 In contrast to conventional quadratic Hamil- 
tonians, though, the number of fermions (not to be con- 
fused with the physical electrons in the system) is not 
conserved, excluding the application of Brouwer's for- 
mula in its existing form. To generalize Brouwer's result 
to this somewhat unconventional case, we follow a path 
similar to the one taken by Vavilov et alZ in studying 
the photovoltaic effect in open chaotic cavities. Starting 
from the nonequilibrium Keldysh Green function tech- 
nique, we show how the adiabatic limit is obtained from 
a systematic gradient expansion. In this manner we are 
able to express the instantaneous spin current in terms of 
an energy-shift matrix^ leading to a Brouwer-type for- 
mula for the adiabatically pumped spin current. 

The formalism outlined above has three notable ad- 
vantages over the scattering approach 24 originally used 
by Brouwer to derive his result: (i) It conveniently ac- 
commodates the case where particles are not conserved; 

(ii) All orders of perturbation theory are summed up in 
the Keldysh technique, thus exceeding linear response; 

(iii) Based on a systematic gradient expansion, one can 
easily read off the small parameter controlling the adia- 
batic limit. We emphasize, however, that the resulting 
Brouwer-type formula for the electronic spin current is 
formally expressed in terms of the scattering matrix for 
the Majorana fermions that appear in the transformed 
Hamiltonian. While technically useful, these degrees of 
freedom have neither a simple representation nor inter- 
pretation in terms of the physical electrons in the leads, 
thus obscuring a clear physical picture. It remains to be 
seen whether a similar expression can be written down 
for the spin current directly in terms of the scattering 
properties of the lead electrons which carry the current. 

The remainder of the paper is organized as follows. In 
Sec. [II] we briefly review the Toulouse limit, introducing 
the different Green functions that will later be used in the 
course of the calculation. In Sec. IIIII we present general 
symmetry considerations and apply them to the problem 
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at hand. In particular, we show that the instantaneous 
charge current is strictly zero in the absence of poten- 
tial scattering, whereas the spin current is zero unless 
the dot couples asymmetrically to the two leads. Pro- 
ceeding with quantitative calculations, we combine the 
Keldysh technique with a gradient expansion in Sec. IIVI 
to derive a Brouwer-type formula for the adiabatically 
pumped spin current in the Toulouse limit. Using this 
formula, a specific class of pumping cycles is analyzed in 
detail in Sec. |V] In particular, we demonstrate a pump- 
ing cycle for which the total spin pumped per cycle is 
closely equal to h, thus operating as a quantized spin 
pump. Finally, we present our conclusions in Sec. IVII 



II. PHYSICAL MODEL AND TOULOUSE LIMIT 

We begin with a brief review of the Toulouse limit, 
and with introducing the different Green functions that 
will later be used in calculating the pumped spin cur- 
rent. The physical system under consideration is shown 
schematically in Fig.[TJ A spin-^ local moment f is em- 
bedded between two leads of noninteracting spin-^ elec- 
trons, undergoing a spin-exchange interaction with the 
local conduction-electron degrees of freedom on either 
side of the junction. As emphasized in the introduction, 
the impurity moment f can either represent an ultrasmall 
quantum dot with a single unpaired electron^ or an ac- 
tual magnetic impurity as in single-atom^ and single- 
molecule^ transistors. 

Since scattering off the impurity is restricted to the 
s-wave channel, one can reduce the conduction-electron 
degrees of freedom that couple to the impurity to one- 
dimensional fields tp aa (x), where a — R, L labels the 
lead (right or left) and a =1,1 specifies the spin orien- 
tations. In terms of the one-dimensional fields, coupling 
to the impurity takes place via the local spin densities at 

the origin: s a p = jE^j'^Lf )^^'^! )' Tnc most 
general form of a spin-exchange Hamiltonian is therefore 

/oo 
ipl a (x)d x ip acr (x)dx 

a=_u,rt. 0=1,4. -°° 

+ E E jf(t)T X S X af3 -» Bgi H(t)T*,(l) 
a,f3—L,R X—x,y^z 

where we have allowed for different exchange couplings 



J" 13 — J\ a , and for a local magnetic field H acting on the 
impurity spin. Here \xb and gi are the Bohr magneton 
and impurity Lande g factor, respectively. Throughout 
the paper we use units for which h = ks = 1, while the 
electronic charge is taken to be — e. Proper units will 
be reinstated in some of the final expressions presented 
below. 

The Hamiltonian of Eq. JT]) is written for general time- 
dependent exchange couplings J" (t) and local magnetic 
field H(t). Our interest, however, will be in slow periodic 
modulations of the transverse couplings J"^(t) = Jy^(t) 
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and the local magnetic field. The longitudinal couplings 
J" 13 will be taken to be constant in time and equal to 
particular values as detailed below. It is this fine tun- 
ing of J"' 3 that defines the Toulouse manifold and which 
enables our exact solution. 



A. Toulouse limit 

The spin-exchange Hamiltonian of Eq. ([T]) is conven- 
tionally derived from the more basic Anderson impu- 
rity model via the Schrieffer- Wolff transformation^ The 
couplings J"' 3 generated in this case are weak, isotropic 
(i.e., independent of A), and satisfy J LL J RR = (J LR ) 2 . 
The Toulouse limit corresponds to a different sector in 
the parameter space of the Kondo Hamiltonian where 

L = J RR = 2ttv f and R = 0. The transverse cou- 
plings Jx^it) = Jy 13 ^) = Jj^{t) and the local magnetic 
field H(t) are allowed to be arbitrary, and will subse- 
quently be taken to be periodically modulated in time. 
Physically, this choice of parameters implies that tun- 
neling is always accompanied by a spin flip. Although 
quite remote from the situation encountered in real quan- 
tum dots, this model is expected to correctly describe the 
strong-coupling regime of the Kondo effect, as argued in 
the introduction and elaborated on in Refs. Ill and 
l37l . In particular, it has been showrt^ that the strong- 
coupling physics of the Anderson impurity model is best 
described both in and out of equilibrium by couplings 
that satisfy 



T LL T RR 

•J 1 J \ 



(J 



LR\2 



(2) 



As described in detail in Ref. [33], the Hamiltonian of 
Eq. (TTJ) can be mapped under the conditions listed above 
onto a free-fermion form. The mapping involves a se- 
quence of steps, comprised of (i) bosonizing the fermion 
fields, (ii) a nonlocal canonical transformation involving 
the conduction-electron spin degrees of freedom, and (iii) 
refermionization of the boson fields to form four new 
fermion fields: ip v (x) with v — c,s,f,sf. Here c, s,/, 
and sf stand for charge, spin, flavor (left minus right), 
and spin-flavor fields. In addition, the impurity spin f, 
which has been mixed by the canonical transformation 
with the conduction-electron spin degrees of freedom, 
is represented in terms of two real Majorana fermions: 
a = —\[2t v and b = —\/2t x . At the conclusion of these 
steps one arrives at a quadratic Hamiltonian conveniently 
written in the form 

= E E e ^l,k^v,k + i^BgiH(t) ba 

iy—c 7 s,f.sf k 

+ iJ?f{t) xtfb + iJ sf{t) XsfU + iJj{t) Xja, (3) 
where we have introduced the three couplings 

ji L {t) + j RR \t) 



(4) 
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Jj{t) = 



j£ L (t)-J™(t) 



(5) 
(6) 



Here the energies ek are equal to —vpk, a is an ultraviolet 
momentum cutoff corresponding to a lattice spacing, and 
L is the effective size of the leads (i.e., k is discretized 
in units of 2ir/L). The fields {v = /, s/) are local 
Majorana fermions, defined as 



- Ipu,) 



iV2L - 



(7) 
(8) 



Relaxation of each of the conditions — 0, J^ L + 

jRR _ ^ Vp ^ anc [ jii _ jflfl _ o introduces a differ- 
ent interaction term into the Hamiltonian of Eq. ^ , as 
discussed in Ref. Hl| and detailed below. 

Although noninteracting, the Hamiltonian of Eq. Q is 
unconventional in the sense that it does not conserve the 
number of ip fermions (not to be confused with the phys- 
ical electrons in the system). Indeed, the fermion fields 
ip v (x) with v = c, s, /, sf have neither a simple represen- 
tation nor a simple interpretation in terms of the origi- 
nal electronic degrees of freedom. Consequently, not all 
observables can be computed based on the mapping of 
Eq. ([3]). Only observables that have a simple representa- 
tion in terms of the ip fields are accessible. Fortunately, 
both the charge and spin currents fall in this category. 

To derive the transformed forms of the electronic 
charge and spin currents, it is necessary to go back to 
their original representation in terms of the physical elec- 
trons in the leads. Denoting the total number operator 
for electrons with spin projection a in lead a by N aiT , the 
charge current flowing from right to left is given by 



I, = 



iH > N L1 + N Ll ] = ie[H,N 



R\ 



N 



Hi J 



(9) 



Here Tt is the Kondo Hamiltonian of Eq. JT]) . Since charge 
fluctuations are excluded on the dot, the instantaneous 
charge current outgoing from the left lead (left commuta- 
tor) is identical to the instantaneous charge current flow- 
ing into the right lead (right commutator). This is no 
longer the case with the spin current, defined as half the 
difference in particle currents between the spin-up and 
spin-down electrons. (The factor of one-half comes from 
the electronic spin projection in the z direction). Indeed, 
the spin currents associated with the left and right leads 
differ by a term proportional to dr z /dt, which stems from 
conservation of the total spin projection S? ota j of the en- 
tire system. Fortunately, this difference in currents has 
no significance for our purposes, since dr z /dt averages to 
zero over a single pumping cycle. This grants us the free- 
dom to work with our operator of choice. In the following 
we shall concentrate on the symmetrized spin current, 
i.e., the average of the spin currents to the left and to 



the right of the impurity, which turns out to be the most 
convenient current combination to work with. With this 
convention, the (symmetrized) spin current flowing from 
left to right is written as 



-[H,N m -N Rl -N M +N Ll ]. 



(10) 



Equations ((9|) and (fT0|) specify the electronic charge 
and spin currents in terms of the physical electrons. The 
transformed operators, I' c and I' s , are obtained by repeat- 
ing the same sequence of steps as applied to the Hamil- 
tonian, namely, bosonization, a nonlocal canonical trans- 
formation, and refcrmionization. Skipping the details of 
the algebra^ we quote here only the end result: 



I'c = ieJ f (*)X/«i 



and 



It 



J sf (t) xt*a- 



(ii) 



(12) 



Note that although these expressions are written in terms 
of Majorana fermions, they describe the actual electronic 
charge and spin currents flowing in the system. The un- 
conventional forms of the currents stem from the nonlocal 
transformation that has been applied. 



B. Keldysh Green functions 

To compute the spin current, we shall make use of the 
nonequilibrium Keldysh Green function technique. The 
basic ingredients of the theory are the greater, lesser, re- 
tarded, and advanced Majorana Green functions, defined 



($(t')a(t)), 



(13) 



(14) 



G££(*.0 =Tie(±tTt')({a(t)J(t')}). (15) 

Here a, (3 € {a, b}, while the upper and lower signs in 
Eq. (fT5| correspond to the retarded (r) and advanced 
(a) Green functions, respectively. The curly brackets in 
Eq. (p~5]) denote the anticommutator. 

In thermal equilibrium, the Majorana Green functions 
are easily found by summing all orders of the perturba- 
tion theory in the time- independent couplings J^j, J J , 
and H. Specifically, switching over to the energy domain 
and assuming the wide-band limit one obtains 



G r ' a {e) 



1 



{t±iT a ){e±iT b )-^ Bgi HY 



iHB9iH e ± iT a 



(16) 
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Here we have adopted a 2 x 2 matrix notation, with the 
indices 1 and 2 corresponding to a and 6, respectively. 
Equation (fT6| features two new energy scales, 



{Jjf + {J; f f 



and 



(17) 



(18) 



where po = 1/(2ttvf) is the density of states per unit 
length in the leads. These two scales determine the 
widths of the various Majorana spectral functions, and 
thus play the role of Kondo temperatures at the Toulouse 
limit. The conventional single-channel Kondo effect is 
best described by the case where r a = Tb = Tjc, corre- 
sponding to the condition specified in Eq. ([2]). The equi- 
librium greater and lesser Green functions are related in 
turn to G r ' a (e) through standard identities: 

i[l-f(e)][G r a0 (e)-G a af3 (e)}, (19) 



G^(e)=i/(e)[G; /3 (e)-G^(e)], 



(20) 



where /(e) is the Fermi-Dirac distribution function. 

As emphasized above, Eqs. (fl6|) — (f20|) are restricted to 
thermal equilibrium. They do not apply when any of the 
couplings jfji J J , and H is time dependent, which is 
the case of interest here. Indeed, time-dependent cou- 
plings are generally difficult to treat analytically even for 
nonintcracting systems. Below we shall first derive the 
instantaneous spin current for a general time-dependent 
setting, but will eventually be interested in slow peri- 
odic modulations of the four coupling constants listed 
above. In terms of the original spin-exchange Hamilto- 
nian of Eq. ([1]), we allow for general time variation of 
the couplings J"' 3 and field H, but demand that the lon- 
gitudinal exchange couplings J"^ be held fixed at their 
Toulouse-limit values. We exclude variations in the phase 
of J± R = (J RL ) , as this corresponds to biasing the 
system. Accordingly, we take J^_ R = J RL to be real 
throughout the paper. 



C. Deviations from the Toulouse limit 

We conclude this section by briefly describing the mod- 
ifications that are introduced into the Hamiltonian and 
the current operators upon departure from the Toulouse 
manifold. As discussed in Ref. [4l|, the Hamiltonian of 
Eq. ^ is supplemented by three new interaction terms 
away from the Toulouse limit: 

H'^H' + H int (21) 



with 



Hi, 



,k' ■ 



k,k' 



(22) 



Here xt are the local Majorana fields of Eqs. ([7]) and 
(|8"|), while : ip 1 fc i/v,fc' : stands for normal ordering with 
respect to the unperturbed Fermi sea of the ip fermions. 
The three couplings J~ = (J^ L - J R R )/2, J+ = (J z iL + 
J RR )/2 — 2-kTlvf and J^ R measure the deviations from 
the Toulouse manifold in each of the three possible direc- 
tions in parameter space. The new tunneling term J^ R 
also modifies the current operators I' c and I' s , which take 
the general forms 



I' c = ieJj(t)xja - ejL R (t)xjxtfba 



(23) 



and 



I' 



^—X f X sf ba. 

(24) 

Here we have explicitly allowed for time variation of the 
new coupling constant J^ R - 



III. SYMMETRY CONSIDERATIONS 

Before proceeding to detailed calculations, in this sec- 
tion we first present general symmetry considerations 
applicable to any two-lead system. By analyzing their 
implications for the Kondo Hamiltonian of Eq. (fT]), we 
identify necessary conditions for finite charge and spin 
currents to be pumped through the system. 



Particle-hole symmetry acting separately on 
each lead 



Consider a general two-lead system where each lead 
is represented by a single spinful channel. The charge 
current flowing into lead a (a = L, R) is given by 



I c ,a = -ie[H,N a1 +N al ], 



(25) 



while the symmetrized spin current I s flowing from left 
to right is specified in Eq. (TTrJ|) . Here N aa denotes the 
total number operator for electrons with spin projection 
a on lead a. Let us consider the situation where the time- 
dependent Hamiltonian, TL, is invariant under a particle- 
hole transformation that converts particles on each lead 
to opposite-spin holes on the same lead (i.e., k a — > 
e «¥w Ca _ fc _^ where a is the spin index opposite to a; the 
phases (p acr are arbitrary) . The total number operator for 
electrons on lead a, N a = +iV Q j, is converted under 
such a transformation to n a — N a , where n a marks the 
total number of electronic states in lead a. Consequently, 
I c _ a transforms according to 



-ie[H, N Q 



-ie[H, n c 



(26) 



k,k' 



If the system begins its evolution from equilibrium, i.e., 
the statistical averaging at time t depends solely on the 
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Hamiltonian at previous times, then the instantaneous 
charge current I c ,a(t) = (Ic.a(t)) = —(Ic,a(t)) must nec- 
essarily be zero. 

The above argumentation is quite general, making no 
reference to the microscopic details of 7i, nor to the tem- 
perature T. Its usefulness lies in revealing the necessary 
(but not sufficient) condition for a finite instantaneous 
charge current to flow: Either the Hamiltonian is not 
permanently invariant under the particle-hole transfor- 
mation indicated above, or the statistical averaging is 
not determined by the Hamiltonian alone (as is the case 
for a finite voltage bias). Note that this symmetry bears 
no information on the spin current, as the latter is in- 
variant under the particle-hole transformation specified 
above. 



B. Particle-hole symmetry that interchanges the 
two leads 

An equivalent statement can be made about the sym- 
metrized spin current I s (t) = (I s (t)) in case of a particle- 
hole symmetry that simultaneously interchanges the two 
leads. Indeed, let us now assume that H. is invariant 
under a transformation where particles on each lead are 
converted to opposite-spin holes on the opposite lead (i.e., 
ot h <7 ^ e lip °"'Ca,-k,s 1 where a is the lead index oppo- 
site to a). Under such a transformation I c . a is converted 
to —I c ,ai while I s is transformed to — I s . Hence, the in- 
stantaneous spin current I s (t) must necessarily be zero 
whenever evolution begins from thermal equilibrium. By 
contrast, no general statement can be made about the 
charge current in this case, apart from the obvious iden- 
tity I c ,a(t) = -I c ,a(t). 

C. Application to the Kondo Hamiltonian 

Our discussion thus far was quite general. We now 
apply the symmetry arguments presented above to the 
Kondo Hamiltonian of Eq. (TIJ. It is easy to verify that 
Eq. (TT]) is invariant under the particle-hole transforma- 
tion 

^U( x ) -> Tpal( x )> W a {( x ) -^a.\{x) (27) 

(corresponding to ijr a k a — > ±ip a ._k,o-), regardless of the 

local field H and the Kondo couplings J?^ = jf Q - 
Hence, the instantaneous charge current for tunneling 
through a Kondo impurity is strictly zero in the absence 
of a voltage bias, as follows from the general discussion of 
subsection IIII A[ In particular, no charge can be pumped 
through the system unless a finite amplitude for potential 
scattering is introduced into the Hamiltonian. Although 
the description of real quantum dots typically requires 
the inclusion of a potential-scattering term, the latter can 
be made negligibly small by operating the device deep in 



the Kondo regime. In this manner charge transport can 
be excluded. 

Similarly, it is straightforward to confirm that the 
Hamiltonian of Eq. |T]) is invariant under the combined 
transformation 

V^t^) ~* V'du(z), i/}' ai (x) -t -ipoer(x) (28) 

(corresponding to ^ aka — ► ±ipa,-k.s), provided the 
intra-lead exchange couplings obey J^ L — J\ R - Thus, 
the instantaneous spin current is strictly zero if the spin 
couples equally to the two leads, as follows from the gen- 
eral discussion of sub-section IIII Bl Spin pumping there- 
fore requires asymmetric coupling to the two leads at 
least in some stretches of time. 

It is instructive to re-derive these results based on the 
symmetries of the transformed Hamiltonian TC + 7ii„t , 
which serves primarily as a check for the correctness of 
Eqs. and (f2"2")l . Other than the free kinetic-energy 
term, the flavor field ijjf enters both Tt' and Hint only in 
the form of x J > which is invariant under the particle-hole 
transformation 

Note that the latter transformation is restricted to the 
flavor sector. Consequently, TH! +7imt is invariant under 
the transformation of Eq. (j2"!)|) , while the charge-current 
operator, being proportional to Xf , transforms according 

to I' c — > — I' c [see Eq. (|23p]. This in turn demands that 
I c (t) be zero in the absence of a voltage bias, in agreement 
with the general symmetry considerations of Eq. (|27p. 

Similarly, when J^ L — J RR the couplings and 
J~ drop from the transformed Hamiltonian H' + 7ii n t , 
which now depends on the field ip s f either through the 
free kinetic-energy term, or in the form of xtf ■ a re " 
suit the transformed Hamiltonian is invariant under the 
spin-flavor particle-hole transformation 

4 f , k -» i>*f,-k, (30) 

while the spin-current operator, being proportional to 
x7/> acc l u i res an extra minus sign: I' s — > — I' s [see 
Eq. ([24]) ]. This in turn implies that the instantaneous 
spin current I s (t) is strictly zero if the leads couple 
equally to the spin, in agreement with the symmetry 
considerations of Eq. (|28p. Interestingly, I s (t) remains 
zero for symmetric coupling also in the presence a finite 
bias, as the latter couples solely to the flavor field. This 
result, originally derived in Ref. [33j for nonequilibrium 
steady state, extends also to time-dependent couplings 
and time-dependent bias. 

IV. PUMPED SPIN CURRENT 

Having established that the instantaneous charge cur- 
rent vanishes for a generic Kondo model in the absence 
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of a voltage bias, we focus our attention hereafter on the 
spin current. To this end, we evaluate I s (t) — (I s (t)} 
exactly on the Toulouse manifold by summing all orders 
of the perturbation theory in the couplings J J (t) , J^j (t) , 
and H{t). We show that a finite spin current can indeed 
be pumped through the system by applying a nonzero 
magnetic field, provided the spin couples asymmetrically 
to the two leads. The calculation proceeds in three steps. 
Using the Keldysh technique, we first derive a formal ex- 
pression for the instantaneous spin current in terms of 
the Majorana Green functions of Eqs. (p~5|) — (p~5]) . This 
portion of the derivation makes no assumption on the 
time-dependent couplings, apart from the restriction to 
the Toulouse manifold and the exclusion of an applied 
voltage bias. The resulting expression is recast in turn in 
subsection IIVBI in terms of a time-dependent scattering 
matrix for the Majorana fermions, defined in Eq. (j4l)|) 
below. The latter scattering matrix reduces in equilib- 
rium to the Fourier transform (with respect to energy) 
of the conventional single-particle scattering matrix. Fi- 
nally, a systematic gradient expansion is carried out in 
subsection II V CI for the case of slowly varying potentials, 
resulting in a Brouwer-type formula for the adiabatically 
pumped spin current. 



A. General formulation 

We begin by formally deriving the instantaneous spin 
current I s (t) using the Keldysh technique, for general 
time-dependent couplings on the Toulouse manifold. As 
is always the case with the Keldysh approach, we as- 
sume that the perturbations J^, JJ , and H have been 
switched on at some distant time in the past, to — ► — oo, 
prior to which the system was in thermal equilibrium. 

To set the stage for the Keldysh formalism, the spin 
current I s (t) is first written as 



where 



j; f (t)G< sf+ (t,t) - j+(t)G< sf _(t,t) 

G< S/+ (M') - (xf f (t')a(t)), 
Gf, sf .(t,t') = (£}(f)S(t)> 



(31) 

(32) 
(33) 



[see Eq. (|12l) ]. Using standard diagrammatics, each of 
the latter correlators is expressed in an exact manner as 



G b,sf- 



where 



+G r ab (t,T)g< <+i+ (T,t')\dT, (34) 

/OO 
■T/W [G< a {t,T)g a sf ^{T,t') 
-oo 

+Gl a (t,r)g< < _ < _(T,t')}dT, (35) 



9?, P ,At,t') = (X P J (t')xl(t)) 



(36) 



and 



gz, p ,At>t')=iW-mm),xt(t')})o (37) 

are the unperturbed Green functions for the local Ma- 
jorana fields. Here v = f,sf and p,p' = ±1. The 
zero subscripts in Eqs. (|36|) and (137)) come to indicate 
that both the time evolution and statistical averaging 
are taken with respect to the unperturbed Hamiltonian, 
i.e., the free kinetic-energy part of Eq. ([3]). 

In writing Eq. (|34p and (f3"S"]) . we have used the fact 

that <7y |_ and gj y identically vanish as long as no 

voltage bias is applied. 43 - Indeed, in the wide-band limit 
Eqs. (|3"6")l and (|3"T|) take the explicit forms 

9%j (*. t') = 27rpoV^(* - (38) 



and 



g^ p , p ,(t,t') = inp 6 pp ,5(t-t'), (39) 

where T(t) is the Fourier transform of the Fermi function 
/(e): 



T(t) = lim / 
»;-o+ J_ { 



2^ 



/W. (40) 



GZ.f+(fJ) = / J tf(-r)[G< b (t,r)g^ +t+ (r,t') 



The limiting procedure used in Eq. (|40j) corresponds to 
regularizing the conduction-electron density of states per 
unit length according to p(e) — poe' 6 '' 1 , and taking the 
wide-band limit D = I/77 — > 00. Equation |39|) is slightly 
modified for a finite bandwidth but remains propor- 
tional to 8 pp i. Inserting Eqs. (J3SJ) and iJSSJ into Eqs. (f34|) 
and (I35p . and plugging the resulting expressions into 
Eq. (|3"l"T) , one obtains 



hit) = i^-J+ f (t)J- f (t) [G<(t, t) - G< (t, t)] + np I [J; f (t)G r ab (t, t) J+ (t) - J+{t)Gl a {t, r)j; f {r)]T{r - t)dr. 

( 41 ) 

It is easy to see at this point that the instantaneous spin current vanishes in the absence of an applied magnetic 
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field, as it physically should. Indeed, setting H — in Eq. Q, the two Majorana fermions a and b decouple within 
the Hamiltonian TC . As a result the Green functions G a b and G ba identically vanish, as does I s . It is also apparent 
that I s is strictly zero unless the impurity couples asymmetrically to the two leads, in accordance with the general 
symmetry arguments of Sec. IIII Cl In fact, I s (t) vanishes not only when J^ L — J± but also for J^ L = — J RR , which 
stems from yet another symmetry of the Toulouse-limit Hamiltonian. Specifically, Eq. ([3]) is invariant for = 

under the particle-hole transformation ip s f,k —> ~Y s f -k> w hile I s transforms according to I s — > —I s . Consequently 
I s {t) = —Is(t) must necessarily vanish when J± = — J RR . 



B. Time-dependent scattering matrix 



Although formally exact, Eq. (|4ip requires knowledge of the time-dependent Green functions G a b and Gf, a , which 
are difficult to compute for a general time-dependent setting. In order to implement the adiabatic limit, it is useful to 
first recast Eq. (|41j) in terms of a time-dependent scattering matrix to be defined below. This goal requires a sequence 
of steps, starting with expressing the lesser Green functions G^ b and G^ a in terms of the retarded and advanced Green 
functions. Since the Hamiltonian of Eq. ([3]) is quadratic, one has the identities 



G<(t,t) = 2-k Pq J™ drJ X dr' G r aa {t, r)| J~ f {r) J~ f (r') + Jj{r) J/ (t')|g^(t', t) 



+ G r ab (t,r)J+(T)J+(T')G a bb (r\t) 



Hr-r'), 



G<(t,t) = 2npoJ °dTj OC dT' Gl a (t,T)^j; f (T)j; f (T') + jy(T)jy(T')^G a aa (r',t) 



+ Gl b (t,T)j: f (T)j: f (T>)G a ba (r',t) 



Hr-r') 



(42) 



(43) 



Substituting Eqs. (|4*2")l and ([4*3")) into Eq. (|4"Tj) . it is convenient to introduce the scattering T- matrix associated with 
the Majorana fields xfj and \J , 



T r > a (t,t') = 2np 



Jtf(t)G r bb a (t,t')Jt f {t') J+(t)G r ba a (t,t>)J; f (t>) J+(t)G r ba a (t 7 t>)Jj(t>) 

J; f (t)G^{t: tVtfW) J^Gl-ait, W^V) J7 f W^(t, t')JJ («') 

jy(t)G r ab a (t,t')jt f (.t') jf(t)G:>2(t,t>)j- f (t') jj{t)G^{t,t<)jj{t>) 



(44) 



Here the row and column indices i = 1,2,3 are identified with (s/, +), (s/, — ), and (/,—), respectively. In terms of 
the T-matrix specified above, the spin current is written as 



/OO pOO pOO 

dr{T r (t, t)T{t - t) - T{t - r)T a (r, *)}+*/ dr dr' T r (t, t)T{t - r')T a (r', t) 
-OO J— OO J —oo 

/OO pOO pOO 

dr{T r {t, t)T{t - t) - T(t - r)T a (r, t)} + i dr dr' T r (t, t)T{t - r')T a (r', t) 
-OO J —OO J —OO 



(sf-,sf+) 



(45) 



(«/+,«/-) 



Finally, the time-dependent scattering matrix for the Majorana fields is defined as 

S(t,t') = S{t-t')l-iT r (t,t r ), 
S^t') = S(t-t')l+iT a (t,t'), 

allowing us to compactly rewrite Eq. (|45[) in the form 



J s (t) = -Im/ dr I dr' S(t,r)F(T - T')S j (r',t') 



(«/+,«/-) 



(46) 
(47) 



(48) 



A word is in order at this point about the time- dependent scattering matrix of Eq. (|46p . Physically, 
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S(t, if) describes the scattering of an incoming Majorana 
fcrmion at time t' to an outgoing Majorana fermion at 
time t. It reduces in equilibrium to the Fourier trans- 
form (with respect to energy) of the conventional single- 
particle scattering matrix, and remains an exclusive func- 
tion of the time difference At = t — t' under general 
steady-state conditions. Although this ceases to be the 
case in the presence of time- varying fields, S(t,t') contin- 
ues to satisfy the generalized unitarity relation 




drS{t,T)S f (r,t') = S{t-t')l, (49) 



to be utilized below. 



C. Gradient expansion and Brouwer-type formula 

The main achievement of Eq. (|4"5)) is the expression 
of the instantaneous spin current in terms of the time- 
dependent scattering matrix S(t, if). For a general peri- 
odic modulation of the couplings Jji, J J , and H, the 
instantaneous spin current at time t depends on the 
specifics of the pumping cycle. For example, the history 



and rates at which parameters are varied. This is not 
the case in the adiabatic limit, where the only informa- 
tion needed to predict the pumped spin per cycle is (i) 
the shape of the pumping trajectory in parameter space, 
and (ii) the equilibrium 5-matrix along the trajectory. 
Similar to adiabatic quantum pumping in noninteract- 
ing systems, the adiabatic limit is approached when the 
characteristic modulation frequency obeys <C r a ,Fb 
at each point along the pumping trajectory. Here T a and 
Tft are the energy scales defined in Eqs. (|17p and l|18|). 
respectively. 

To substantiate these claims and devise a Brouwer- 
type formula for adiabatic quantum spin pumping in the 
Kondo regime, we resort to a systematic gradient ex- 
pansion of Eq. (|48|) . To this end, we first introduce the 
Wigner transform of the time-dependent scattering ma- 
trix, 

S(e,T)= £V e - T s(T+^r-0. (50) 

Next we apply the well-developed machinery of the Gra- 
dient expansion^ For example, the Wigner transform of 
the convolution of two functions, 



dre leT i dr'A [T 

-oo J — OO 



')B (t',T- 



(51) 



has the formally exact representatio: 



,45 



[A * B](e, T) = e ^{ d Td?-d?df}A(e, T)B{e, T) = A(e, T)B(e, T) + i (d T Ad e B - d e A0 T B) + ■■■. (52) 

Here d A and d B stand for differential operators that act on A(e,T) and B(e,T), respectively. The usefulness of 
Eq. ([52]) comes into play when the expansion on the right-hand side is controlled by a small parameter. This is 
indeed the case in the present context, where the double convolution of Eq. (j48|) possesses an analogous expansion in 
gradients of S(e, T). Each combined derivative drd e is parametrically reduced for Eq. (|4"5)) by a factor of il/T, where 
r is some characteristic value of either r a or Tb in the relevant time interval. The scale T is bounded from below by 
the minimum of T a and Tf, along the pumping cycle, a quantity denoted hereafter by T. Hence, for f2 <C T one can 
settle with linear order in dxd e to obtain 



I,(t) =Im 



47T 



SSt + - [(ftS)(0,St) - (0 e S)(ftSt)] \ f(e) + - {S(9 t St) - (ftS)S+} 



2/ 



df(e) 



(«/+.«/-) 



r 



(53) 



All terms omitted in this expression are of order (fl/T) 2 
or higher, and thus can be safely neglected4£ 

The term proportional to the Fermi function /(e) in 
Eq. (fB"3")) is purely diagonal to order 0(SI/T), as can be 
seen by expanding the unitarity relation of Eq. (|49|) to 
first order in time gradients: 

SSt + I [($S)(0 e St) - (d e S)(d t &)] +0[(Q/T) 2 ] = 1. 

(54) 

Since Eq. (|53|) requires an off-diagonal matrix element of 



the expression in the square brackets, the instantaneous 
spin current reduces in the adiabatic limit to 



d( 



I s (t) = Re —f'(e) 

, 07T 



s(ftst) - (a t s)st 

L J («/+,»/-) 

This expression can further be simplified by noting that 
S(e,t) is equal to leading order in Q/T to the instanta- 
neous scattering matrix, i.e., the equilibrium scattering 
matrix with all system parameters J s * , J J , and H frozen 
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at their instantaneous values at time t: 



S(e,t)=S cq (e;J±(t),JJ(t),H(t))+0(n/r). (56) 



Consequently, one can substitute S oq in for S in Eq. (J55J). 
Lastly, one can exploit the unitarity of the equilib- 
rium S'-matrix, S oq Sj q = 1, to replace (<9 t S cq )Sj q with 
— S cq (9 t S < t q ) in Eq. (|55j) . This yields the final expression 
for the spin current, 



I s {t) = Re 



4tt 



/'W [S(^St)] 



(«/+,»/-) 



(57) 



Here and in the remainder of the paper the symbol S 
is used as a shorthand for the instantaneous scattering 
matrix S oq (e; J± (t), Jj(t), H(t)), 

Equation (|5T[) is exact in the adiabatic limit, £1 — > 
0. Its derivation was based on a systematic truncation 
of higher order terms in Q, controlled by the expansion 
parameter Q/T. It therefore encompasses all pumping 
trajectories and all coupling regimes, whether weak or 
strong. This should be contrasted with the commonly 
used linear-response theory, which is restricted, strictly 
speaking, to weak coupling only. 

In the following we shall consider examples of pumping 
cycles where two system parameters, generically termed 
X\ and X2, are varied slowly and periodically in time 
along a certain closed trajectory C in parameter space. 
The quantity of interest in this case is the total magne- 
tization in the z direction, or spin, transferred from left 
to right in a single pumping cycle. The latter quantity is 
defined as 



(S) = J>I s (t)dt, 



(58) 



where I s (t) is the instantaneous spin current. Using 
Eq. (f57|) one can express (S) as a line integral along the 
contour C, 



<«> = 



(«/+,«/-) ' d — 



This expression applies to the variation of any number 
of system parameters X 1; • • • , Xjy. In the particular case 
where N = 2, one can make use of Green's theorem to ex- 
press the spin pumped per cycle as a geometric property 
of the Majorana-fermion scattering matrix, analogous to 
Brouwer's formula for noninteracting systems. Explicitly, 
(S) assumes the form 



(S) = Re J —/'(e) J^dX 1 dX 2 



(«/+,»/-) 



,(60) 



where A is the (oriented) area in parameter space en- 
closed by the contour C. 

Equation (|60|) is the central result of our study. We 
devote the remainder of the paper to analyzing its im- 
plications for a particular class of pumping trajectories 
defined below. 



x 2 

h 



Y(b) 



: (c) 
-->■ 



<- 

(a) 



\Jsf 

(d)A 



-h 



FIG. 2: The pumping cycle under consideration in Sec. [V] 
The first pumping parameter, Xi, controls the Kondo cou- 
plings J~j and J J , which vary according to J~* — X\ and 



J 



The third Kondo coupling, u 3 p 



held fixed throughout the cycle, along with Tk- The sec- 
ond pumping parameter, X2, controls the Zeeman splitting 
fisgtH, which varies according to fiBgiH — Xi. 



V. APPLICATIONS 

We conclude our analysis by applying the formula 
derived above to study a particularly simple class of 
pumping cycles where one parameter, X\, controls the 
transverse Kondo couplings, and the other parameter, 
X2, controls the applied magnetic field. Based on our 
Toulouse-limit calculations we will show that such a cy- 
cle can be used to realize a pure quantized spin pump, 
namely, quantized spin pumping without any charge 
transport. 

To make contact with realistic systems such as quan- 
tum dots, we impose hereafter the condition J^ L jf R — 



(J± R ) 2 , corresponding to (J s f) 2 



{Jf ) — _ \ u sf) 

As mentioned above, this condition best describes the 
strong-coupling physics of the Anderson impurity model, 
where a single Kondo scale T a = Tb = Tk emerges. 
Keeping J^, and thus Tk, fixed, we parameterize 
JJ , and HB9iH according to 



J 



and 



HBgiH = X 2 . 



(61) 
(62) 

(63) 



In terms of the original Kondo couplings to the two leads, 
Eqs. (foTj) and (|6"2"|) translate to 



LR 



xi 



(64) 
(65) 
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The pumping cycle under consideration is depicted 
schematically in Fig. [51 It consists of four segments, two 
in which X\ is tuned from ±Ji to ^fJjf while X 2 is 
kept fixed [lines (a) and (c)], and two in which X 2 is 
tuned from ±h to ^fh while Xi is held fixed [lines (b) 
and (d)]. The cycle C thus consists of periodic open- 
ing/closing of the transverse couplings to the left/right 
leads, followed by inversion of the applied magnetic field 
at points where spin-flip scattering is restricted to one 
lead only. The analogous cycle for real quantum dots 
comprises of periodic opening/closing of the tunneling 
rates to the left/right leads, followed by inversion of the 
applied magnetic field at points where tunneling is re- 
stricted to one lead only. 

Combining Eq. §Q for (S) with Eqs. (|gg ]l . (|g7 ]) . (P J) . 
and (|16p for the instantaneous S'-matrix, one obtains af- 
ter some straightforward but tedious algebra 



7T 



def'(e) 



Re 



hT, 



K 



V + (T K - ief 
h + e 



arctan 



(66) 



Here we have restored Ti for proper units of (S). Repre- 
sentative plots of (S) as a function of both h and T are 
shown in Fig. [3] As expected of the Kondo regime, (S) 
is an exclusive function of the rescaled parameters h/Tx 
and T /Tk- In particular, at T = one finds 



(S) 



2Ti 

IT 



T k h 



arctan 



h 

Tk~ 



. (67) 



which has the formal expansion (S) /% = l—O [(Tk /h) 3 ] . 
Hence, the pumped spin per cycle is closely quantized to 
h when the magnetic field H performs a large enough ex- 
cursion along the pumping cycle. The effect of a temper- 
ature is to reduce the spin pumped per cycle. However, 
(S) remains closely quantized to h when h 3> Tk , T. Im- 
portantly, when T <C Tk, it suffices that h will only mod- 
erately exceed Tk in order for (S) to closely approach %. 
For example, at T = the pumped spin per cycle is equal 
to 0.82ft (0.967i) by the time h — Tk (h — 2T K )- 

The above results were derived at the Toulouse limit, 
which does not correspond to any realistic parameters. It 
is therefore pertinent to question the relevance of these 
results to actual quantum dots. Since any exact solution 
can be used to extract universal low-energy properties of 
the Kondo effect, we expect Eq. (|66|) to be quantitatively 
correct when T, h <C Tk ■ Equation (|6"6"]) should remain 
qualitatively correct as one of the parameters T or h be- 
comes comparable to Tk, though quantitative deviations 
are expected. Still, since (S) approaches h quite rapidly 
with increasing h (essentially by h ~ Tk), and since the 
departure from strong coupling is only logarithmically 
slow in h, we expect (S) to remain nearly quantized in 
real quantum dots provided T <C Tk- This picture is 
further supported by a naive application of Brouwer's 
formula using the exact T = single-particle scatter- 
ing matrix^i and by slave-boson mean-field theory of 
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FIG. 3: The spin pumped per cycle, (S), in units of %. (a) 
Plotted as a function of H/Tk for different T; (b) Plotted as 
a function of T '/Tk for different h. For T < Tk, the spin 
pumped per cycle rapidly approaches h with increasing h. 
Explicitly, (S) exceeds 0.9R for all T < 0.452V when h = 2T K . 



the corresponding Anderson models The Toulouse limit 
fails, however, to describe the weak-coupling regime, as 
certain bare couplings are required to be large. In partic- 
ular, Eq. (I66p should neither be quantitatively nor qual- 
itatively correct when Tk <C T. 

A simple interpretation of Eq. (|66p follows from the 
observation that the ground state of the Kondo model 
is that of a local Fermi liquid. Only resonant elastic 
scattering takes place at the Fermi level when T = 0, 
as reflected in the Abrikosov-Shul resonance. The latter 
resonance is pinned to the Fermi energy when H = 0, 
and is split by an applied magnetic field. This basic phe- 
nomenology can be mimicked by a simple noninteracting 
resonant-level model, 



Hrlm = ^2 ^ek^l^kaa ~ (J-BgiH(d\d^ -d\di) 

a=L,R k,cr 

+E ^{v>La+ h - c -}> (68) 



which is studied below. Here i\)\ aa creates an electron 
with wave number k and spin projection a on lead a 
(a = L,R), while d\ creates a localized electron on the 
level. 

Allowing for slow periodic modulation of H and V a in 
Eq. (|68p . we extract the adiabatically pumped spin and 
charge along a closed pumping cycle analogous to the 
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one shown in Fig. [2] For a generic trajectory C in the 
parameter space (Xi,X 2 ) defined below, the adiabati- 
cally pumped spin and charge are given for Hrlm by the 
standard Brouwer formula 

— /'(e) J dX 1( iX 2 
x [d Xl S a d X2 Si - d X2 S a d Xl St] LL (69) 

and 

^J_ 00 2tti J a 
x [d Xl S a d X2 St - d X2 S a d Xl St] LL . (70) 

Here a =f,4 and a = ±1 are used interchangeably to 
label the spin projection. The domain of integration, 
A, is the (oriented) area in parameter space enclosed by 
the contour C. The instantaneous S'-matrix pertaining to 
Wrlm is written in the L-R basis as 



S ff (e) = 



where 



1 - 2iT L Gl{e) -2 l VrZr^G^(e) 
-2iVTir^(e) l-2ir fl G£(e) 



1 



e - onBgiH + ir H 



is the associated dot Green function. Here T 4 



(71) 



(72) 



-r, 



with r Q = np V 2 is the resonance width, which plays the 
role of the Kondo temperature in the Kondo model. 

By analogy with the cycle of Fig. [21 we vary the two 
pumping parameters X\ = Tl — Tr and X 2 = psgiH 
while r + is held fixed. As before, the cycle is composed 
of four segments, two in which X\ is tuned from ±r + 
to =Fr+ while X 2 is kept fixed, and two in which X 2 
is tuned from ±h to =ph while X\ is held fixed. Using 
Eqs. (|69"|) - (|72)) for this cycle one obtains (Q) = and 



7T 



def(e) 



Re 



hT, 



h 2 + (r+ - it) 2 
h + e 



arctan 



(73) 



Both results are identical to those obtain at the Toulouse 
limit, provided r + is identified with Tk ■ Thus, the phys- 
ical picture underlying Eq. ([55)) is consistent with that 
of simple resonant elastic scattering, where a single res- 
onance is symmetrically split about the Fermi energy by 
an applied magnetic field. 

From a theoretical standpoint it is clear that one can 
realize a quantized spin pump using either a quantum 
dot in the Kondo regime or a Zeeman-split single-particle 
resonance that is tuned to the Fermi energy. However, 
practical considerations make the Kondo-dot scenario a 
more promising candidate for the realization of such a 
device. Indeed, modulation of the couplings to the two 



leads is typically accompanied in real devices by a capac- 
itive shift of the dot level. In case of a simple resonance, 
the induced modulation of the dot level will generally 
produce a finite charge current, and is likely to spoil the 
quantization of the pumped spin. The Kondo-dot sce- 
nario is immune to such fluctuations, as these produce 
only a tiny shift of the Abrikosov-Shul resonance. In- 
deed, as discussed in Sec. IIII1 charge transport is strictly 
forbidden as long as the Coulomb-blockaded dot can be 
described in terms of a pure Kondo Hamiltonian having 
no potential scattering. Although a realistic description 
of quantum dots generally requires the inclusion of po- 
tential scattering, the latter term can be made negligibly 
small by operating the device deep in the Kondo regime. 
In this manner charge transport can be excluded. 



VI. CONCLUSIONS 

In this paper we have presented an exact analysis of 
adiabatic quantum pumping through a quantum dot in 
the Kondo regime. It follows from general symmetry ar- 
guments that the instantaneous charge current is strictly 
zero in the absence of potential scattering and for zero 
voltage bias. A similar statement applies to the sym- 
metrized spin current either in the absence of an applied 
magnetic field or for symmetric coupling to the leads. 
Pumping of a spin current therefore requires both a fi- 
nite magnetic field and for left-right symmetry to be si- 
multaneously broken. Both conditions are readily met 
in practical devices, making ultrasmall quantum dots a 
natural candidate for the realization of a spin battery. 

To quantify this statement, we have computed the 
pumped spin current exactly at the Toulouse limit. Ex- 
ploiting the mapping onto a quadratic Hamiltonian and 
performing a controlled expansion in the small parameter 
fl/Tx (fi being the characteristic modulation frequency, 
Tk is the Kondo temperature), we have expressed the 
pumped spin per cycle as a geometric property of the 
scattering matrix associated with three flavors of Majo- 
rana fermions, which are the effective degrees of freedom 
at the Toulouse limit. In particular, employing the cou- 
pling to the leads as one pumping parameter and the 
applied magnetic field as another, we have shown that 
one can devise pumping cycles that realize a pure quan- 
tized spin pump. Namely, a device for which the average 
spin pumped per cycle is closely equal to h, but where no 
accompanying charge current is produced. We expect the 
pumped spin per cycle to remain nearly quantized in real 
quantum dots provided that one operates at T <C Tk ■ 

There have been by now a number of different pro- 
posals in the literature for the realization of spin pumps, 
employing diverse setups such as chaoticity in quantum 
dots?i£ ferromagnetic leadsf^ spin-orbit interactions ) 16 ! 19 
classical turnstile cycles^ one-dimensional Luttinger- 
liquid physics ^± and finally the Kondo effect in quan- 
tum dotsJ^ While all these proposals reported schemes 
to realize a pure adiabatic spin pump along specific cy- 
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cles, the quantization of the spin pumped per cycle has 
been shown to be the case only in the classical turn- 
stile setup 1 ^ and for a Luttinger liquid.— In contrast to 
Ref. Il8l the pumping scheme investigated in this paper 
offers an interesting possibility to realize a coherent quan- 
tized spin pump, in which the absence of charge current is 
essentially warranted along all possible cycles (including 
beyond the adiabatic limit). 

The quantization of the pumped spin per cycle re- 
ported in this paper is subject to small deviations as the 
temperature T becomes of order Tjc, or as the magnetic- 
field excursion is altered. Moreover, it applies only to 
the average spin pumped per cycle. In order to better 
characterize such a quantum pump, a detailed study of 
its noise properties (and full counting statistics) is desir- 
able. A study of the statistical properties of the Kondo 



pump is a challenge left for future work. 
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